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ABSTRACT 


Context. 

Aims. We perform an extensive characterization of the broadband emission of Mrk 421, as well as its temporal evolution, during the non-flaring 
(low) state. The high brightness and nearby location (z=0.031) of Mrk 421 make it an excellent laboratory to study blazar emission. The goal is to 
learn about the physical processes responsible for the typical emission of Mrk 421, which might also be extended to other blazars that are located 
farther away and hence are more difficult to study. 

Methods. We performed a 4.5-month multi-instrument campaign on Mrk 421 between January 2009 and June 2009, which included VLBA, F- 
GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. This extensive radio 
to very-high-energy (VHE; £ > 100 GeV) y-ray dataset provides excellent temporal and energy coverage, which allows detailed studies of the 
evolution of the broadband spectral energy distribution. 

Results. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves 
show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) 
at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise 
behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes 
with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the 
non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. 
Conclusions. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi¬ 
wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring 
during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk 421 during non¬ 
flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within 
the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton 
scenario, indicates that the electron energy distribution is most variable at the highest energies. 

Key words. Galaxies: BL Lacertae objects: individual: Markarian 421 
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1. Introduction 


Blazars are a class of radio-loud active galactic nuclei (AGN) 
where the relativistic jet is believed to be closely aligned to our 
line of sight. They emit radiation over a broad energy range from 
radio to very high energy y rays (VHE; £ > 100 GeV), which 
is highly variable at all wavelengths. Their spectral energy dis¬ 
tributions (SED) are dominated by the jet emission and show 
two bumps, one at low energies (radio, optical. X-rays) and the 
other at high energies (X-rays, y rays, VHE). While the origin of 
the low-energy bump is presumably synchrotron emission from 
relativistic electrons, the origin of the high-energy bump is still 
under debate. To constrain current theoretical models for broad¬ 
band blazar emission, simultaneous observations of those ob¬ 
jects over the whole wavelength range and over a long period 
are needed. It is important to perform observations at typicaQor 
even lower states in order to have a baseline to which other (flar¬ 
ing) states can be compared, as distinct physical processes might 
play a role when the source is flaring. Weak blazars in a low state 
are particularly poorly studied in y rays because of the difficulty 
to detect them at these energies with current instrumentation. In 
addition, most multi-wavelength programs are triggered when a 
source is flaring, and not when it is in low state. 

The high-energy peaked BE Lac object (HBL) Mrk421 was 
the first extragalactic object discovered at VHE ( [Punch et al.| 
1992| l. It is one of the brightest extragalactic X-ray/VHE ob¬ 
jects, and because of its proximity (z = 0.031) the absorption 
by the extragalactic background light (EBL) is low ( [Albert et al.| 
2007[ l. Mrk421 has been well-studied during phases of high ac¬ 
tivity, but simultaneous broadband observations in a low state, 
covering both energy bumps, were missing until recently. 


Starting in 2009, a multi-wavelength (from radio to VHE), 
multi-instrument was organized to monitor the broadband emis¬ 
sion of Mrk421. The scientific goal was to collect a complete, 
unbiased and simultaneous multi-wavelength dataset to test cur¬ 
rent theoretical models of broadband blazar emission. In this pa¬ 
per we analyze the temporal variability of Mrk421 in all wave¬ 
lengths during the 4.5-month observation period in 2009. During 
the entire period, Mrk421 did not show any major flaring activ¬ 
ity (e. g., [Aleksic et al.[[2012[ Fortson et al.[[2012| [Fossati et al. 


2008 Gaidos et al.|1996[ [Mankuzhiyil et ai.[[2011 1. The multi¬ 

wavelength dataset is used to enhance our understanding of the 
origin of the high-energy emission of blazars beyond the usually 
observed flaring states. The underlying physical mechanisms re¬ 
sponsible for the acceleration of particles in jets are compared 
with those observed during flares. This paper can be understood 
as a sequel to Abdo et al. (201 lb I ( [Paper I[ l, where the SED of 
the Mrk421 2009 data was analyzed. 


This paper is organized as follows; Sect.|^introduces the par¬ 
ticipating instmments and the multi-wavelength data. The anal¬ 
ysis of the variability in each waveband is presented in Sect. 
Cross-correlations and periodic behavior are examined in Sects. 
1^ and 1^ and Anally in Sect. we summarize and discuss our 
results. 


* The complete data set shown in Fig.[^is available in electronic form 
at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or 
via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 

' We use the term “typical” instead of “quiescent”, to describe a 
state that is neither exceptionally high/flaring, nor at the lowest possible 
level. Even though the term “quiescent” has been used in the past to de¬ 
note non-flaring activity in Mrk421 and other blazars, we note that the 
term quiescent refers to the lowest possible emission, which is actually 
unknown, and hence not suitable in this context. 


2. The 2009 multi-wavelength campaign 


The duration of the 2009 campaign on Mrk421 was 4.5 
months from 2009 January 19 (MID 54850) to 2009 June 1 
(MJD 54983). 29 instruments participated in the campaign. The 
intended sampling was one observation per instrument every two 
days, whenever weather, technical and observational limitations 
allowecQ The list of participating instruments and the time cov¬ 
erage as a function of energy range are shown in Table 2 and 
in Figure 6 of [Paper I The schedule of the observations can be 
found onlin^H The individual datasets and the data reduction are 


presented in detail in section 5 of Paper I[ and will therefore not 
be introduced again in this paper. Besides the datasets reported 
in Paper I, this paper also reports VHE data from 115 hours 
of dedicated Mrk421 observations with the Whipple 10-meter 
telescope (operated by the VERITAS collaboration). These data 
are essential for the excellent temporal coverage in the VHE for 
this campaign. Details on the light curve presented here can be 


found in Pichel (2009) with the general Whipple analysis tech¬ 
nique described in [Horan et al.[ ( [2007[ l and [Acciari et al.[ ( [2014j l. 
The frequencies/wavelengths covered by the campaign are ra¬ 
dio (2.6 - 225 GHz), near-infrared (7, H and K), optical (B, V, 
g, R and /), UV (Swift/UWOT Wl, W2 and M2), X-ray (0.3- 
195 keV), high-energy (HE) y rays (0.1-400 GeV) and VHE 
(0.08 - 5.0 TeV). 

Results on the broadband SED as well as a detailed discus¬ 


sion of the SED modeling can be found in Paper I It is the most 
detailed SED collected simultaneously for Mrk421 during its 
typical activity state and the first time where the high-energy 
component is completely covered by simultaneous observations 
from the Fenni-LP^T and the VHE instrument MAGIC. This al¬ 
lowed the characterization of the typical SED of Mrk421 with 
unprecedented detail. In Paper I the SED could be modeled rea¬ 
sonably well using either a one-zone synchrotron self-Compton 
(SSC) model having two breaks in the electron spectrum, or a 
hadronic (synchrotron proton blazar, SPB) model. In order to 
distinguish between these two scenarios, one must look at the 
multi-wavelength variability. One- and multizone SSC models 
predict a positive correlation between X-ray and VHE flux vari¬ 
ations (e.g., Graff et al.|2008[), as they are produced by the same 


electron population. In the SPB models of Paper I a strict cor¬ 
relation between those two bands is neither generally expected 
nor excluded, but can appear when electrons and protons are 
accelerated together. Furthermore, the one-zone SSC model of 
Paper I[ predicts a correlation between low-energy y rays from 
Fermi-LAT with millimeter (from SMA) and optical frequen¬ 
cies, something which would be hard to incorporate in the SPB 
model, as the radiation is produced at different sites. 

In the following sections we will first characterize the flux 
variability in all wavebands and then have a detailed look at 
the cross-correlation functions between light curves of different 
bands, primarily at X-rays vs. VHE and optical vs. X-rays, HE 
and VHE correlations, but also at all other combinations as they 
might reveal something interesting. 


^ e.g., for imaging air Cherenkov telescopes (lACTs) like MAGIC 
or Whipple, observations during moonlight are only possible to a very 
limited extent, resulting in regular gaps of ~ 10 days in the VHE light 
curves. 

^ https://confluence.slac.stanford.edu/display/ 
GLAMCOG/Campaign+on+Mrk421+(Jan+2OO9+to+May+2809) 
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Fig. 1. Light Curves of Mrk421 from radio to VHE from 2009 January 19 (MJD 54850) to 2009 June 1st (MJD 54983). Vertical 
bars denote flux measurement errors, and the horizontal bars denote the time bin widths into which some of the light curves arg 
binned. The Ferm/-LAT photon fluxes are integrated over a three-day-long time interval. The Whipple 10-meter data (with an energy 
threshold of 400 GeV) were converted into fluxes above 300 GeV using a power-law spectrum with index of 2.5. 
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log(frequency [Hz]) 


Fig. 2. Fractional variability F^ar as a function of frequency. 
Open circles denote F^ar values in f?-band calculated with the 
host galaxy subtracted as prescribed in|Nilsson et al.|p007|l 


3. Variability 

3.1. Light curves 

Figure[T]shows the light curves from radio to VHE. No substan¬ 
tial (larger than a factor of 2) flaring activity happened during the 
campaign; however, some level of variability is present in all en¬ 
ergy bands. In the radio band the variability is least pronounced. 
A significant level of variability is present in the near-infrared 
(NIR), optical and UV accompanied by an overall increase in 
flux with time. At X-ray, HE and VHE there is also consid¬ 
erable variability, and only a small overall downward trend in 
the overall X-ray and VHE flux with time is observed. The X- 
ray flux variations are stronger on average than the variations 
in the other wavebands, but still much weaker than the maxi¬ 
mum values historically registered for the X-ray and VHE bands 
( [Balokovic et al.|2013| [Cortina & Holder|2013| l. 


3.2. Fractional variability 


In order to quantify and characterize the variability at different 
energy bands, we calculated the fractional variability 


A 


(o-L) 




( 1 ) 


i.e., the excess variance normalized by the flux, according to 
Vaughan et al. ( 2003| l, where S is the standard deviation of the 
N flux measurements, is the mean squared error and (x)^ 

is the square of the average photon flux. We estimate the uncer¬ 
tainty of Fvar according to|Poutanen et al. (20081, 


•y ^ vai 


• err(o-2 ). 


( 2 ) 


where err(cr^„„) is given by equation 11 in Vaughan et al. 

( |2^ : - 


err(cr: 


Nxs) “ 


A 


2 \C^e, 


N <x>2 



(3) 


This prescription to calculate the uncertainties is more appropri¬ 
ate than equation B2 in [Vaughan et al.| (|2003|l for light curves 


that have an error in the excess variance comparable to or larger 
than the excess variance. This is, however, not the case for most 
light curves in our sample, as AFyar according to Poutanen et al. 
(|2008|l is less than 5% smaller compared to equation B2 in 


Vaughan et al.| ( 200^ . Eor the Fenni-LALL, and SwiftfBALT light 
curves, the difference is 10%. 

The Fvar values and errors for the different energy bands (in¬ 
struments) are shown in Eigure]^ As already noticed when look¬ 
ing at the light curves, Mrk 421 shows little variability in radio, 
and low but significant variability in all other wavebands with 
the largest variability in X-rays. We note that in the 2-10 keV 
band it is intrinsically more variable than in the 0.3-2 keV band, 
a characteristic which has been recently reported for Mrk 421 
during high X-ray and VHE activity ( Aleksic et al.pOlSbj l 

Because of the instrument sensitivity, the Ferm/-LAT, 
SwiftlBALT and RXTE/ASM data are binned into 3- and 7-day 
bins (instead of 1-day bins) and therefore, the variability on 
smaller timescales is not probed, so Fy^ might be underesti¬ 
mated. When rebinning the RXTE/PCA light curve (sampled 
every ~ 2 days) into 7-day bins, Fvar decreases by ~ 15% 
and agrees with Fvar for RXTE/ASM within the errors. The 
Swift/XRT light curves are irregularly sampled. There were mea¬ 
surements every ~ 7 days during the early part of the campaign 
but there are also large gaps and a period of sub-daily observa¬ 
tions. Rebinning the SwiftPCRT light curves into 1- and 7-day 
bins does not change Fy^y by more than a few percent and all 
values agree within the errors. 

The results reported in Eigurej^are not affected by the tem¬ 
poral binning or the uneven sampling of the light curves and 
hence can be considered as characteristic of Mrk 421 during the 
multi-wavelength 2009 campaign (see Appendix A for details). 

It is interesting to compare these results with the ones re¬ 
ported recently for Mrk 501 in Doert et al. (|2013|l and Aleksic 


et al. (2015a|l, where the fractional variability increases with 


energy and is largest at VHE, instead of X-rays. The compar¬ 
ison of these observations indicates that there are fundamental 
differences in the underlying particle populations, environment, 
and/or processes producing the broadband radiation in these 
two archetypical VHE blazars. The higher X-ray variability in 
Mrk 421 might also be related to the higher synchrotron domi¬ 
nance with respect to the one observed in Mrk 501. According 
to the broadband SEDs measured for Mrk 501 and Mrk 421 
during the typical (non-flaring) activity (Abdo et al. 201 la|b| l, 

4 X 


^^SynCpeak ^ 


2 X for Mrk 501 and vFf > 


^^iCpeak ]y[j.]j 42 i These SEDs were parametrized within the 
one-zone SSC framework in [Abdo et'af] ( |201 la|b| l, using, for 
Mrk 421, a magnetic held B ~ 2.5 times higher than the one 
used for Mrk501 (38 mG vs. 15 mG), which naturally produces 
a synchrotron bump that is relatively higher than the inverse- 
Compton bump. The higher magnetic held in Mrk 421 may also 
lead to a higher variability in the X-ray band (with respect to the 
y-ray bump) through a faster synchrotron cooling of the high- 
energy electrons (Xcooi-Sync V-B^)- 


3.3. Evolution of the X-ray spectral shape with the X-ray flux 

Systematic variations of the X-ray spectral shape are a com¬ 
mon phenomenon during blazar flares (e.g. Fossati et al.|2000] l. 
A harder-when-brighter behavior is quite typical during flares 
in blazars, and this characteristic has already been observed in 
Mrk 421 (e.g. Tramacere et al.|2009 1 . Sometimes one can iden¬ 
tify loops in the photon index vs. flux diagram during the course 
of a flare, which could be related to the dynamics of the sys- 


4 














































































J. Aleksic et al.: The 2009 multiwavelength campaign on Mrk 421: Variahility and correlation studies 




Fig. 3. Left: X-ray hardness ratio for the SwiftfXKI bands 2-10 keV and 0.3-2 keV vs. the X-ray flux in the 2-10 keV band. Right: 
Power-law index of RXTE/PCA spectra above 3 keV vs. the X-ray flux in the 2-10 keV band. 


tem, as reported by Kirk & Mastichiadis ( 1999| l or Rieger et al. 
( 2000| l. Such behavior was also observed in Mrk 421 during a 
big flare in 1994 ( |Takahashi et al.|ri996[ ). Here we investigate 
whether these patterns exist in Mrk 421 during its typical (non¬ 
flaring) activity. 

The SwiftpiKY spectrum cannot be fit with a simple power- 
law because this instrument covers the peak of the synchrotron 
bump, and hence the X-ray spectrum in the 0.3-10 keV band is 
curved. We can quantify the hardness of the 5w//f/XRT spec¬ 
tra by using the ratio of the X-ray fluxes in the bands 2-10 keV 
and 0.3-2 keV, and study its evolution with respect to the X-ray 
flux in the 2-10 keV band. This is shown in the left-hand panel 
of Fig. The RXTE/PCA spectra, which cover the falling seg¬ 
ment of the synchrotron bump (when the source is not flaring), 
can be fit with a simple power-law function, and hence here we 
can report the spectral slope vs. the X-ray flux in the 2-10 keV 
band. This is shown in the in the right-hand panel of Fig. Both 
Swift and RXTE show clearly that the X-ray spectra harden when 
the X-ray emission increases, and hence we can confirm that the 
harder-when-brighter behavior also occurs when the source is 
not flaring. We also investigated the temporal evolution of the 
plots shown in Fig. looking for loop patterns in the spectral 
shape-flux plots (clockwise or counter-clockwise) but we did not 
And any. 


3.4. Power density spectrum 

Another way to characterize the variability of a given source 
is the power spectral density (PSD). The PSD quantifies the 
variability amplitude as a function of Fourier frequency (or 
timescale) of the variations. The derivation of the PSD is based 
on the discrete Fourier transform of the light curve under con¬ 
sideration. For blazars, the shape of the PSD is usually a power- 
law Py cx with spectral index a between 1 and 2 (Abdo 
et al.|201^ Chatterjee et al.|2012| l, i.e., there is larger variabil¬ 
ity at smaller frequencies / longer timescales. This is generally 
referred to as “red noise”. Other features in the PSD, such as 
breaks or peaks, indicate characteristic timescales or (quasi)- 
periodic signals. 

Calculating the PSD via the discrete Fourier transform is 
straightforward for light curves that are frequently and regularly 
sampled over a long period of time. However, in reality, observa¬ 
tion time is usually limited and often interrupted by bad weather. 


object visibility and technical issues, i.e., we are normally deal¬ 
ing with unevenly sampled light curves of limited length that 
may have large gaps, and this has serious effects on the measured 
PSD. If the light curve is discretely sampled instead of contin¬ 
uous (which is usually the case as we are dealing with discrete 
observations or values that are binned over a certain time pe¬ 
riod), its Fourier transform is convolved with a windowing func¬ 
tion, which becomes very complicated when a light curve has an 
uneven sampling and gaps ( |Merrifleld & McHardy|199^ . In ad¬ 
dition, light curves of a finite length are affected by red-noise 
leak, i.e., variability below the smallest frequency (or largest 
timescale) probed. This variability power leaks into the observed 
frequency band and changes the observed PSD shape. This effect 
can manifest as a rise/fall trend throughout the entire time inter¬ 
val of the light curve. Likewise, aliasing, i.e., variability power 
from frequencies larger than the Nyquist frequency, affects the 
variability in the observed frequency range. The effect of sam¬ 
pling on the study of periodicities will be discussed in Section 

m 


These effects of the sampling pattern on the PSD can be 
avoided by using the simulation-based approach of|Uttley et al. 


( |20021 l (PSRESP) to derive the intrinsic PSD of a light curve 
and its associated uncertainties. We applied this Monte Carlo fit¬ 
ting technique following the prescription given in the appendix 
of Chatterjee et al. ( |2008) l to all light curves with ~ 30 or more 
flux measurements. 

First we generated a large set of simulated light curves using 
the method of |Timmer & Koenig] ( |1995] l. In order to accommo¬ 
date the problems introduced by the sampling of the light curve, 
the simulated light curves were about 100 times longer than the 
measured light curve and then clipped to the required length. 
This way, they suffer from red-noise leak in the same way as 
the measured light curve. The simulated light curves also had 
a much finer sampling than the measured light curve and were 
then binned to the required binning to include the aliasing effect. 
Finally, the simulated light curves were resampled with the ob¬ 
served sampling function, so that the windowing function is the 
same for the measured and the simulated light curves. Poisson 
noise was added to each simulated light curve to account for 
observational noise. As a model for the underlying PSD of the 
simulated light curves we assumed a power-law shape and varied 
the power-law index a in the range 1.0 to 2.5 in steps of 0.1. We 
generated 1000 simulated light curves per a value and per mea- 


5 

































J. Aleksic et al.: The 2009 multiwavelength campaign on Mrk 421: Variahility and correlation studies 


sured light curve. We then calculated the PSD for each measured 
light curve and the 1000 simulated light curves of each model, 
taking as the PSD the modulus squared of the mean subtracted 
light curves’ discrete Fourier transform between the minimum 
frequency Vmin = IIT and the Nyquist frequency VNy - NjlT. T 
is the duration of the light curve. The frequency range covered 
by our data is approximately 10^^- > 10“^ ^ s“^ (corresponding 
to ai 1/120 days^^ - « 1/6 days^'), differing somewhat between 
the light curves depending on the individual length and binning. 
The light curves were binned into 2 - 7-day bins, depending on 
the light curve characteristics. The goodness-of-ht of each PSD 
model w as determined according to the recipe given in the ap¬ 
pendix of|Chatterjee et al. |2008|: The observed function 


2 (PSDobs - PSD,i^) 

(APSD,i™)2 


(4) 


from the observed PSDobs, the average of the 1000 PSDs 
from simulated light curves PSDsin), and the standard deviation 
APSDsim was compared to the simulated distribution 


2 

Tdist. 


>inax 

St., = £ 


(PSDsim,/ - PSDsim)^ 

(APSD,i„)2 


(5) 


calculated from each of the 1000 individual PSDs from simu¬ 
lated light curves PSDsim,;- The success fraction (SuF) is then 
the fraction of^^ larger than ^2^^. 

Figure [4| shows the SuFs for selected measured light curves 
and Table gives the best-ht power-law indices a and their 
uncertainties, calculated as the half width at half maximum 
(HWHM) of the SuF vs. a curve. 

We tried a range of binnings and found that the location 
of the maximum does not vary systematically with the binning 
(Figure]^. The uncertainties, however, depend strongly on the 
light curve binning and on the logarithmic binning of the PSD 
( jPapadakis & Lawrence|1993| l. We used light curve bin sizes of 
a few days (between 2 and 6) and a factor of 1.2 or 1.3 by which 
the logarithmically spaced frequency bins are separated, in order 
to reduce the scatter in the PSD points. Sometimes the flux mea¬ 
surements have large uncertainties (mean error non-negligible 
compared to the variance of the light curve as, e.g., in the case 
of Fermi-LP^) or the light curve has large gaps and/or relatively 
few data points (e.g., MAGIC). In these cases a is mostly uncon¬ 
strained. Large gaps or very uneven binning may result in large 
SuF differences (e.g., MAGIC), or even in changes of the over¬ 
all shape (e.g., OVRO). In these cases there is no good way to 
bin the data and obtain a reliable a. If the binning is too small, 
large gaps are hlled with interpolated (probably unrealistic) data 
in order to calculate the discrete Fourier transform. This results 
in unwanted changes in the reconstructed PSD and thus in unre¬ 
liable a values. If the binning is so large that it accomodates also 
the large gaps such that not too much interpolation is necessary, 
there are too few data points left and the covered frequency inter¬ 
val becomes too small for a reasonable PSD ht. Unfortunately, 
the light curves are too short to make analyses based on contigu¬ 
ous parts of a light curve between large gaps. 

The maximum of the SuF vs. a curve is generally > 0.8 for 
all light curves, but we note that for certain binnings the maxi¬ 
mum SuF can be signihcantly lower, or saturate at 1.0. Thus a 
power-law seems to be a reasonably good ht for all light curves, 
but as explained above, the ht is restricted by the limited fre¬ 
quency range and by the light curve sampling or gaps. In Table 
[T]we give the a values where the SuF has a maximum, i.e., the 


Table 1. Power spectral density (PSD) index a and half-width 
at half-maximum of the success fraction (SuF) for light curves 
with more than 30 hux measurements. 


Instrument 


HWHM 

OVRO 

2.0 

1.2-2.3* 

GRTI 

1.5 

0.7 

MITSuME Ic 

1.6 

0.7 

GASP 

1.9 

0.5 

GRTR 

1.4 

0.8 

GRT V 

1.5 

0.5 

MITSuME g 

1.5 

0.6 

GRTB 

1.4 

1.0 

UVOTWl 

1.4 

0.5 

Swift/XKT (0.3-2 keV) 

1.5 

0.6 

SwiftPiKF (2-10 keV) 

1.4 

0.6 

RXTE/PCA 

1.4 

0.6 

Fermi-hAT 


b 

1 

bo 

MAGIC 

1.6 

0.9 

Whipple 

1.3 

0.6 


Notes. 

slope of a power-law P(v) oc 

® The SuF vs. a curve is asymmetric, thus the PSD of OVRO is not 
very well constrained. This might he due to the very low variability in 
radio compared to the other wavelengths. 

The PSD of Fermi-hAT is largely unconstrained because of the 
large flux measurement errors. The SuF is approximately constant and 
high (>0.8) over a large a range with no clear maximum. 


best-htting power-law indices. As uncertainties we mention only 
the median HWHM of the distribution of HWHM from differ¬ 
ent light curve binnings between 2 and 6 days, but please note 
that this value itself has an uncertainty. In many cases it is not 
clear why we should prefer one binning over another, so a cer¬ 
tain range in SuF shapes is possible. It should be pointed out that 
these uncertainties in deriving the width of the SuF vs. a curve 
do not affect the analysis in the following sections, as we always 
use the best-htting a, which does not vary with the binning. 

The best-htting PSD models for most light curves are found 
to be power-laws with indices ~ 1.3 - 1.6. There are no big dif¬ 
ferences between the SuF vs. a curves of different instruments. 
The only exceptions are GASP R band and OVRO, where the 
shape is different and a is higher (though still in agreement with 
the other a values within the uncertainties). A possible explana¬ 
tion might be the small hux error bars and the dense sampling 
(several values per night, i.e., larger frequency range compared 
to all other light curves). Without logarithmic binning and with 
a light curve binning of 1 or 2 days there is a clear maximum at 
1 . 8 . 


There is no evidence of a break in the (relatively short) fre¬ 
quency range covered by the multi-wavelength data. Simulated 
light curves with an underlying broken power-law PSD did not 
improve the success fraction. The X-ray PSD power-law indices 
are similar to the ones reported in Kataoka et al. ( |2001[ ) for the 
same frequency range. The PSD shape and a are consistent with 
what was found for blazars by other authors (|Chatteijee et al. 
20081 X-rays; [Abdoet al.|2010| HE y rays). 


4. Cross-correlations 


We use the discrete cross-correlation function (DCF) method of 


Edelson & Krolik (1988 1 to quantify the correlation of the flux 


variations between all possible light curve pairs, as long as the 
light curve has more than 30 flux measurements, i.e., we use 
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Fig. 4. Success fraction as a function of power-law index a of the PSD for three selected light curves (GRT V band, RXTE/PCA 
and MAGIC) and a range of light curve binnings (2-6 days) to illustrate the effect of the binning. The location of the maximum 
does not change significantly with the binning, but there is significant variation in shape, width and amplitude. 


the same set of light curves as in Section 3.4 This way we can 
assess correlations between VHE, HE y rays. X-ray, UV, optical 
and some radio frequencies. 

As a cross-check, we calculate for each light curve pair also 
the z-transformed cross-correlation function (ZDCF; pklexander] 
|1997| l. The ZDCF is based on the DCF, but the bin widths 
of the ZDCF are chosen such that the number of points is 
the same for each bin, i.e., they are different-sized as opposed 
to the DCF, where all bins have the same size. In addition, 
Fisher’s z-transform is applied to the cross-correlation coeffi¬ 
cients. According to |Larsson| ( 2012| l, the ZDCF is more ro¬ 
bust than DCF for undersampled (w.r.t. the flux variations) light 
curves. However, for well sampled light curves, the ZDCF has 


lOkeV) bands is shown in Figure]^ The correlations peak at time 
lag Af = 0 and appear to be strongly significant (>5cr), when 


been shown to be consistent with the DCF (e.g., Dietrich et al. 
|1998| Smith & Vaughan]|2007| l. For this study we used mostly 


the DCF, as the temporal bin is fixed, so that it also allows us to 
trivially compare and even combine results from different pairs 
of instruments. Therefore, we used the ZDCF for verification 
purposes only. 

Two different approaches are used to determine the uncer¬ 
tainties of the DCF. The easiest and fastest way is to simply use 
the errors given in Edelson & Krolik (1988 1 . However, as dis¬ 
cussed in lUttley et al.| ( |20031 l, these are not appropriate for de¬ 
termining the significance of the DCF when the individual light 
curve data points are correlated red-noise data. Depending on 
the PSD and the sampling pattern, the significance as calculated 
by Edelson & Krolik ( |1988| l might be overestimated. To get a 
better estimate on the real significance of the correlation peaks 
we used a Monte Carlo approach, following the descripion of 


Arevalo et al. (2009 1 . The Monte Carlo technique is described in 
detail in Section 14711 

We used a binning of six days for all DCFs because this way 
different DCFs can be easily compared or, if needed, combined. 
To make sure that we do not miss correlations or time lags, we 
always tried a range of binnings, depending on the sampling of 
the involved light curves. 

As the following paragraphs will show, significant correla¬ 
tions are only found between X-rays and VHE. In addition, X- 
rays and optical light curves seem to follow opposite trends. 


4.1. VHE —X-rays 

The correlation of the flux variations between the VHE 
(MAGIC, Whipple) and X-ray (RXTE/PCA and SwiftfXRT 2- 


considering only the errors calculated according to [Edelson & 
|Krolik|p988| l (black error bars). However, as mentioned above, 
the error bars calculated using Edelson & Krolik ( |1988| l can 
overestimate the real significance of the correlation. To get a bet¬ 
ter estimate on the real significance of the correlation peaks we 
use the following Monte Carlo approach. 

For each X-ray light curve we created a set of 1000 simulated 


light curves in the same way as in Section 3.4 using a power- 


law with the best-fitting slope as determined in Section [374| from 
the PSRESP method (see Table [T]). The X-ray flux was sam¬ 
pled more often and with smaller statistical errors than the VHE 
flux, and hence, in order to ascertain the confidence levels in the 
DCF calculation, it is reasonable to use simulated RXTE/PCA 
and SwiftPiRT (2-10 keV) light curves instead of the VHE 
light curves. We cross-correlated each of the 1000 simulated X- 
ray light curves with the observed MAGIC and Whipple light 
curves. A power-law spectmm with index 2.5 ( |Hillas et al.|1998| l 
was used to normalize the integral flux of the Whipple 10 m 
data (with an energy threshold of ~ 400 GeV) to an integral flux 
above 300 GeV in order to provide a comparison with MAGIC. 
We therefore also cross-correlated the original and the simu¬ 
lated X-ray light curves with the combined Whipple-i-MAGIC 
light curve. From the distribution of 1000 DCFs (i.e., the cross¬ 
correlations of the simulated X-ray light curves with the real 
VHE light curves) we then calculated the 95 and 99% confi¬ 
dence limits, and show them in Figure For each combina¬ 
tion of RXTE/PCA and SwiftfXRT (2-10 keV) with MAGIC, 
Whipple and Whipple-i-MAGIC the DCF shows a peak at time 
lag At - 0 with a probability larger than 99%. There are no other 
peaks or dips in the DCF between VHE and X-rays that appear 
significant. A positive correlation between X-rays and VHE has 
been reported many times during flaring activity (e.g., [Fossati 
|et al. 12008] ), but has never been observed for Mrk 421 in a non¬ 
flaring state. Our simulations show that the real significance of 
the correlation is 3 - 4cr, indeed confirming that the error bars 
calculated using Edelson & Krolik ( |1988| l slightly overestimate 
the significance of the correlation. 

We average the DCFs, which has the advantage that spuri¬ 
ous features are smoothed out while features present in all DCFs 
(i.e., those features that are real) are strengthened. This is partic¬ 
ularly useful when having many possible combinations and/or 
marginally significant features like the ones that will be reported 
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At (days) 



At (days) 



At (days) 


At (days) 


Fig. 5. The DCF of the combined Whipple and MAGIC (VHE) light curve, correlated with the RXTE/PCA and the SwiftPiKT 
(2-10 keV) light curves are shown in the upper left and lower left panels. The black error bars represent the uncertainties as derived 
from |Edelson & Krolik| ( |198^ . The green lines represent the 1% and 99% extremes of the DCE distribution of simulated RXTE/PCA 
light curves when correlated with the measured VHE light curve. The blue lines represent the 5% and 95% extremes. The upper 
right panel shows the average of the VHE-RXTE/PCA and VHE-5w//f/XRT (2-10 keV) DCEs, with the corresponding confidence 
intervals derived from averaging the DCEs of the simulated light curves. See text for details in the calculation of the average DCEs 
and contours. The lower right panel shows the z-transformed DCEs. 


in Section |43] The DCFs were averaged in the following way: 
For a number of ^ -H 1 real light curves A, Bi, ..., Bq we first 
calculate all q correlation functions DCF(ABi), ..., DCF(AB^) 
using a binning of 6 days. Then we calculate the average DCF 

_ 1 9 

DCF = - y DCF(AB,)- (6) 

There is no prescription to combine the uncertainties derived 
from |Edelson & Krolik] ( fT988[ l for several DCFs, thus no error 
bars are shown in the upper right panel of Figure For the de¬ 
termination of the averaged confidence limits we correlate the 
simulated light curves ai, ..., or bq i, ..., bq^„ (n = 1000, 
q e 1,... ,q) with the original light curve Bq or A using either 

_ 1 9 

DCFsin,,, = - y DCF(a,B,) V; = 1,..., n (7) 

or 

_ 1 9 

DCF,i^,,- = - y DCF(Af7,.,,.) V; = 1,..., n. (8) 


From this distribution of n averaged correlation functions 
DCFsimj we then compute the 95 and 99% confidence limits. 
Whether we use equation |7] or equation depends on the sam¬ 
pling of the light curves. The sampling and statistical uncertain¬ 
ties of the X-ray light curves are much better than the sampling 
of the VHE light curves. Therefore the PSD could be better con¬ 
strained in the X-ray case and hence the confidence limits ob¬ 
tained from correlating the original VHE with simulated X-ray 
light curves are more reliable than the confidence limits obtained 
from correlating the original X-ray with simulated VHE light 
curves. Thus here the light curve A is Whipple-i-MAGIC, the 
light curves B,- are RXTE/PCA and Swift/XRT (2-10 keV) and 
we use equation]^ 

The upper right panel of Figure shows the DCF averaged 
over the two combinations Whipple+MAGIC - RXTE/PCA and 
Whipple-i-MAGIC - SwiftjXKI (2-10 keV) and the correspond¬ 
ing confidence limits. There is a clear correlation at time lag 
Af = 0 with a high confidence > 99%. 

The ZDCEs between the VHE and the X-ray light curves 
show the same behavior as the corresponding DCEs. However, 
as the binning is different for each ZDCE (ranging from sub-day 
scales around time lags At ^ 0 days to a few days at time lags 
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At (days) At (days) 


Fig. 6. DCF of the combined Whipple and MAGIC (“VHE”) 
light curve, correlated with the Swift/XRT (0.3-2 keV) light 
curve. The black error bars represent the uncertainties as de¬ 
rived from |Edelson & Krolik ( |1988 1 . The green lines represent 
the 1% and 99% extremes of the DCF distribution of simulated 
Swift/XRT light curves when correlated with the measured VHE 
light curve. The blue lines represent the 5% and 95% extremes. 


Fig. 7. The average DCF of the Fenni-LAT HE y-ray light curve, 
correlated with the optical and UV light curves (GASP /?-band, 
GRT BVRI, MITSuME g, MITSuME Ic and UVOT Wl), is 
shown in black. The green lines represent the 1% and 99% ex¬ 
tremes of the likewise averaged DCE distribution of simulated 
optical/UV light curves when correlated with the real Fenni- 
LAT light curve. The blue lines represent the 5% and 95% ex¬ 
tremes. 


At ss 60 days), it is not possible to combine them as we did with 
the DCFs. Rebinning the ZDCFs to even 6-day bins, averaging 
them subsequently and using the simulated light curves to assess 
the uncertainties yields almost identical results to the averaged 
DCF. 

Both VHE and X-ray light curves show a weak negative 
trend with time. To make sure that this trend is not responsi¬ 
ble for the correlation, we calculated the (z)DCFs also for the 
detrended light curves. The difference is marginal. In addition, 
when comparing X-ray and VHE light curves, one can see that 
the light curve features nicely agree, i.e., the long-term trend has 
only a minor contribution to the correlation peak, which is driven 
by shorter timescale variability. 

Eigure shows the DCE of the combined Whipple and 
MAGIC light curve, correlated with the Swift/XRT (0.3-2 keV) 
light curve. Although the DCE has a peak at time lag Af 0 days 
which seems to be significant with a confidence level of around 
5cr when considering the Edelson & Krolik ( |1988| l errors, the 
simulations show that this level of correlation is not significant 
(only 1.5cr). 


4.2. HE y rays — UV/optical 

Eigure|7]shows the combined DCF of the Fermi-LAT light curve, 
correlated with optical and UV light curves (using simulated op¬ 
tical and UV light curves to estimate the uncertainties|^ There 
is a peak in the DCF at a time lag At - 0 days, but it is not 
significant. The uncombined DCFs also show a small peak with 
a significance <3cr or no peak at all when using the Edelson & 
|Krolik|p988| l uncertainties. The significance of the small peaks 
is even lower (<2cr) when using the uncertainties from simu¬ 
lated data. Larger light curve binsizes would reduce the errors, 
but also lead to significantly fewer data points, a reduced PSD 
frequency range, and increased DCE bin sizes. Given the small 


In Section 


3.4 


we showed that it is not possible to constrain the PSD 


of the Fermi-LAT light curve because of the large error bars of the light 
curve data points. 


time window (4.5 month long campaign) under consideration, 
we cannot improve the DCE quality by rebinning. 

4.3. X-rays — UV/optical 

We also searched for correlations between the X-ray and 
UV/optical bands. As done in the previous subsections, we 
calculated the DCEs for all possible combinations between 
RXTE/PCA and SwiftfXRT with Swift/UWOT, GASP R band, 
GRT BVRI, MITSuME g and MITSuME Ic. The results are 
shown in panel A of Figure One feature that is common in 
almost all DCFs is an anti-correlation at a time lag At -20 
to -10 days, i.e., optical/UV variations lead X-ray variations by 
~ 15 days. This feature is significant above 99% and is con¬ 
firmed by the ZDCF However, it is not immediately clear what 
might cause this anti-correlation. The first thing that becomes 
apparent when looking at the light curves is the long-term trend. 
The UV/optical light curves show a strong positive trend, while 
the X-ray light curves show a slight negative trend. Therefore it 
is not surprising that the DCFs show an overall anti-correlation 
spread over a large range of time lags. However, this characteris¬ 
tic cannot explain the above-mentioned (anti-)correlations with 
time lags of 10-20 days. 

Hence we detrended the light curves by fitting and sub¬ 
tracting a first-order polynomial to each light curve and re¬ 
calculated the DCFs and the ZDCFs. They are shown in panel 
B of Figure in comparison to the correlations of the origi¬ 
nal light curves (panel A). In the detrended light curves we find 
two results: 1) The overall negative correlation spread over most 
time lags disappears. 2) Some peaks become evident at time lags 
At ^ -36, -18, h- 6 and h- 18 days (the latter two being absent in 
the RXTE/PCA - UV/optical DCF). In the ZDCFs these fea¬ 
tures are generally less pronounced. The presence of such fea¬ 
tures leads to the suspicion that an underlying quasi-periodic 
behavior may be responsible for the (anti-)correlations. Indeed 
there are several local peaks and minima in both the X-ray and 
UV/optical light curves. In Figure|^we illustrate how well these 
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features correlate by overplotting two light curves (SwiftPiKT 
(0.3-2 keV) and GASP R band). Both light curves are normal¬ 
ized. The GASP light curve is also rescaled such that both light 
curves cover approximately the same normalized flux range. In 
addition, the GASP light curve is shifted in x direction by -36, 
-18, 6 and 18 days. For time lags where an anti-correlation was 
detected (-18 and h- 18 days), we also flipped the GASP light 
curve about the horizontal axis (such that light curve peaks be¬ 
come troughs and vice versa), such that in each panel of Figure[9] 
we should see that both light curves follow the same path when¬ 
ever there is a real (anti-correlation) present. However, it is ob¬ 
vious from these plots that some, but not all of these features 
are loosely correlated (as indicated by the low statistical signifi¬ 
cance of the DCF peaks) and that the limited time window ham¬ 
pers the ability to detect a convincing correlation. The behavior 
illustrated in Figure]^ may well happen just by chance without 
being caused by an underlying physical mechanism. In Section 
l^we show that there is no hint of a periodic signal in any of the 
light curves. 

4.4. VHE — optical/UV 

The VHE light curves, when correlated with UV and optical light 
curves, produce a strong negative peak at time lag Af = 0 days 
(Figure[T0]l, i.e., they are anti-correlated with a probability larger 
than 99%. However, the optical/UV light curves show a strong 
positive trend while the VHE light curves show a weak negative 
trend. After detrending the light curves, the anti-correlation at 
time lag At - 0 days disappears. Instead, the DCE now shows a 
similar behavior as the X-ray-optical DCFs (marginally signifi¬ 
cant anti-correlation at time-lag At -18 days and correlation 
at -36 days). This is not surprising given the positive corre¬ 
lation between VHE and X-ray fluxes. More data in the typical 
state are needed in order to judge whether this behavior is just a 
chance (anti-)correlation or if it is caused by underlying physical 
mechanisms. 


light curves each. The simulated light curves were produced in 
the same way as in Section 3.4 and have the same underlying 
PSD (estimated above with PSRESP) and the same sampling as 
the original light curve. From the distribution of LSPs we de¬ 
termined the 95% and 99% confidence limits. We did not And 
significant ESP peaks in any of the light curves. A peak around 
P 18 days is present in a few optical and X-ray LSPs, but al¬ 
ways below 99% confidence level, and in most LSPs even below 
95% confidence level. 


5.2. Autocorrelation 

We use the discrete correlation function ( jEdelson & Krolik] 
|1988| ) to calculate the discrete auto-correlation function (DACE) 
of the variability of Mrk 421 in all observed wave bands. Equally 
spaced and repeated features in the DACF might be a hint to 
characteristic timescales and quasi-periodicities. As in the pre¬ 
vious sections, we use simulated light curves to estimate the sig¬ 
nificance of DACE peaks, i.e., the observed light curve is corre¬ 
lated with 1000 simulated light curves. This results in confidence 
limits that are not symmetric around zero, although the DACE it¬ 
self is symmetric. The origin of the asymmetry relies on the pro¬ 
cess used to determine the confidence intervals, which uses 1000 
Monte Carlo realizations of one light curve, together with the 
asymmetry of some of the light curves. This results in different 
Monte Carlo realizations when the light curve is truncated on the 
left or on the right (negative or positive time lags), hence yield¬ 
ing different results for the confidence intervals. Consequently, 
the asymmetry in the confidence intervals is particularly strong 
where the sampling and variability of the light curve changes 
significantly with time (e.g., UVOT). Figure [m shows DACFs 
for a few representative light curves. There are secondary peaks 
in some DACFs. However, they are all well below the 95% limit, 
i.e., they do not appear to be significant. Hence we do not And 
significant periodicities or characteristic timescales. 


4.5. Other correlations 

All the data in the optical and UV bands vary simultaneously as 
shown in Figure The NIR bands seem to be well correlated 
with the optical and UV bands. However, the number of flux 
measurements per NIR light curve is too small to calculate a 
meaningful DCF or ZDCF 

No significant correlations are found between radio or HE y 
rays with other wavelengths. 


5. Periodicities 


5.1. Lomb-Scargle periodogram 

Although the Lomb-Scargle periodogram (ESP) is not the best 
way to determine the PSD for red-noise data ( Kastendieck et ah] 
201 l| l, it is a good way to And periodicities when dealing with 
unevenly sampled light curves. 

A peak in the ESP at a certain time lag can mean that there is 
a periodicity. However, the sampling also produces peaks, e.g., 
if there is a flux measurement every 2 days, there will be a strong 
peak at period P - 2 days. Uneven sampling may result in one 
or more peaks, if at least part of the flux measurements fol¬ 
low an approximately regular observation schedule. The ESPs 
were determined for periods < L/5 a; 25 days, where L is the 
length of the light curve. To estimate the significance of poten¬ 
tial ESP peaks, we also calculated the ESP for 1000 simulated 


6. Discussion of the main observationai resuits 


Even though Mrk 421 is known for extreme X-ray and VHE vari¬ 
ability, with short and intense flares (e.g. Aleksic et al.||201^ 
Fossati et al.|2008[ Gaidos et al.|1996 1, the X-ray and VHE ac¬ 
tivity measured in the 2009 observing campaign was relatively 
mild, with X-ray/VHE flux variations typically smaller than a 
factor of 2. The VHE flux of Mrk 421 was also relatively low, 
with an average flux of about 0.5 times the flux of the Crab 
Nebula, which is typical for this source ( |Acciari et al.||2014) l. 
Regardless of the low activity, Mrk 421 showed significant vari¬ 
ability in the portions of the electromagnetic spectrum where 
it emits most of its energy power, namely optical/UV, X-rays 
and HE/VHE y rays. The optical/X-ray bands bring informa¬ 
tion from the rising/falling segments of the low-energy bump, 
while the HE/VHE bands tell us about the rising/falling seg¬ 
ment of the high-energy bump. As reported in Section 3.2 (see 
Eigurej^, the highest variability occurs at X-rays {F^ar ~ 0.5), 
then VHE {F,.ar ~ 0.3), and then optical/UV/HE (F„„, ~ 0.2). 
It is interesting to compare these results with the ones reported 
recently for Mrk501 ( Aleksic et al.||2015a[ Doert et ar]|2013| l, 
where the fractional variability increases with energy and is 
largest at VHE, instead of X-rays. The comparison of these two 
observations indicates that there are fundamental differences in 
the underlying particle populations and/or processes producing 
the broadband radiation in these two archetypical VHE blazars. 
Within the one-zone synchrotron self-Compton scenario, which 
is commonly used to model the emission of VHE blazars, the 
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A) not detrended: 





B) detrended: 
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Fig. 8. (A) The DCFs of each X-ray light curve (SwiftfKKT (0.3-2 keV), SwiftPiKT (2-10 keV) and RXTE/PCA (2-10 keV)), 
correlated with several optical and UV light curves, are averaged over all optical to UV bands and shown in black in the upper 
panel. The green lines represent the 1% and 99% extremes of the likewise averaged DCF distribution of simulated optical/UV light 
curves when correlated with the observed X-ray light curve. The blue lines represent the 5% and 95% extremes. The lower panel 
shows the z-transformed DCFs, which were, for the purpose of direct comparison with the DCF, rebinned to the same binning as 
the DCFs and averaged in the same way. (B) shows the same as (A), but all light curves have been detrended (as described in |4.3| l 
before correlation. 


X-ray and VHE variability is driven by the dynamics of the pop¬ 
ulation of relativistic electrons through their synchrotron and 
inverse-Compton emission, respectively. Within this scenario, 
and for typical model parameters, the ~keV emission is domi¬ 
nated by higher-energy electrons, whereas the ~100 GeV emis¬ 
sion is produced by a mixture of lower-energy electrons that 
inverse-Compton scatter in Thomson regime, and high-energy 
electrons that inverse-Compton scatter in Klein-Nishina regime 


(see |Pape7T| . Consequently, the multi-band fractional variability 
reported in Section 3.2 indicates that the population of higher- 
energy electrons varies more than that at lower energies. 


It is worth noticing that the fractional variabilities de¬ 
tected in the energy range 2-10 keV measured by RXTEfPCA, 
RXTE/ASM and Swift/XRT agree reasonably well with a value 


11 























J. Aleksic et al.: The 2009 multiwavelength campaign on Mrk 421: Variahility and correlation studies 


1.6 

1.4 

1.2 

1 

0.8 

0.6 

1 . 6 - 

1.4- 

1 . 2 - 

1 - 

0 . 8 - 


At=-36 days 


i I 'i/ 


I %*' 


I 


'i ^ ' 

y s. 

Ill ■ 


I 


54860 54880 54900 54920 54940 54960 54980 55000 

MJD 


I 


- i 


0.6: # 

r , i 


J 


^ At=-i-6 daySj 
i 

i 1 

T^i" 1 “i 

it , ^ I ‘ If .■ 

OiT I i i jl I r I 

‘ 'L'i ^ 

■ ■ ■ .T,- i 


54860 54880 54900 54920 54940 54960 54980 55000 

MJD 



MJD 


X 

u. 


1 . 6 - 

1.4- 

1 . 2 - 

1 - 

0 . 8 - 

0 . 6 - 




ii'i 


At=+18 days 


f 

Til T i S' II I I I 

': ■ ’ 1 *, 


. i*' 


54860 54880 54900 54920 54940 54960 54980 55000 

MJD 


Fig. 9. Example plots to illustrate the correlation between the optical and the X-ray flux variations. All four panels show the nor¬ 
malized Swift/XRT (0.3-2 keV) light curve in black. The GASP /?-band light curve, normalized and rescaled to match the same 
flux range as the SwiftpiRY light curve, is overplotted in red with different time lags Af = -36, -18, -i-6, and +18 days. In case of 
anti-correlation (Af = -18 and +18 days), the GASP light curve is also flipped vertically. 
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Fig. 10. The average of the DCFs of the combined Whipple and MAGIC (VHE) light curve, correlated with several optical and UV 
light curves, are shown in black in the left panel. The green lines represent the 1% and 99% extremes of the likewise averaged DCF 
distribution of simulated optical/UV light curves when correlated with the observed VHE light curve. The blue lines represent the 
5% and 95% extremes. The right panel shows the same as the left panel, but all light curves have been detrended (i.e., fitted and 
subtracted by a first-order polynomial) before correlation. 


of 0.4-0.^ despite the different observing windows of these 
three different instruments. On the other hand, the fractional 
variability measured by SwiftjXKY in the energy range 0.3-2 
keV is ~0.25, which is a factor of 2 lower than the variability 
detected by SwiftpiRY in the 2-10 keV energy range. Given that 
these two observations are performed with the same instrument, 
the difference in the fractional variability cannot be ascribed to 


As discussed in Section 


3.2 


the somewhat lower value of 
RXTEjASM, 0.33+/-0.03, is due to the 7-day integration time, which 
prevents the detection of variahility with temporal scales of days. 


a different observing temporal period that might include or ex¬ 
clude a particular flux variation, and hence the higher variabil¬ 
ity in the 2-10 keV energy range, in comparison to that in the 
0.3-2 keV energy range, is an intrinsic property of Mrk 421 dur¬ 
ing the 2009 observing campaign, which has also been recently 
reported for Mrk 421 during high X-ray and VHE activity ob¬ 
served in 2010 ( [Aleksic et al.|2015bl l. Because the characteristic 
synchrotron frequency of relativistic electrons is proportional to 
the square of the energy of the electrons (Vc E^), the higher 
synchrotron energies will be produced by higher energy elec- 
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Fig. 11. The discrete auto-correlation function for a few light curves are shown in black. The green lines represent the l%and 99% 
extremes of the likewise averaged DACF distribution of simulated light curves when correlated with themselves. The blue lines 
represent the 5% and 95% extremes. 


trons, and hence this further supports the theoretical framework 
of higher variability in the number of higher energy electrons. 
Given that the high energy electrons are the ones losing their en¬ 
ergy fastest {TcooI *). in order to keep the source emitting 
X-rays, injection (acceleration) of electrons up to the highest en¬ 
ergies is needed. We therefore conclude that the injection (accel¬ 
eration) of high-energy electrons is likely to be the origin of the 
flux variations in Mrk 421. 

From the multi-instrument light curves shown in Figure [T] it 
can be seen that, while the variability in X-ray and VHE oc¬ 
curs mostly on ~day timescales, the variability at OTtical/UV 
occurs mostly on ~week or even longer timescale^ The dif¬ 
ferent variability timescales do not show up in the results re¬ 
ported in Section 3.4 (see Table [^. However, this might be the 
result of the limited sensitivity of the PSD analysis due to the 
uneven sampling of the light curves, and the rather small range 
of frequencies sampled 10^^ ^ s'*), which do not pro¬ 

vide a long enough lever arm to determine accurately (and ulti¬ 
mately to distinguish) the index of the power-law spectrum of the 
PSDs from the different energy bands. Therefore, the multi-band 
light curves and fractional variability show similarities in the 
X-ray and VHE flux variations, which differ from characteris¬ 
tics of the optical/UV flux variations. This observation is further 


® Because of the 3-day span and the relatively large statistical un¬ 
certainty in the flux points, we cannot evaluate whether short or long 
timescales dominate the variability in the HE y-ray band. 


confirmed by the cross-correlation results reported in Section ^ 
which show a positive correlation with no time lag between the 
X-ray and VHE emission, but not between the optical/UV and X- 
ray or VHE. This result indicates that both the X-ray and VHE 
emissions are co-spatial and produced by the same population of 
high-energy particles. It is worth noting that, while such a corre¬ 
lation has been reported many times for Mrk 421 during flaring 
activity (e.g. [Fossati et al. 2008 1 , this is the first time that this 
is observed during a non-flaring (typical) state. Therefore, to¬ 
gether with the observed harder-when-brighter behavior in the 
X-ray spectra, we interpret this experimental observation as evi¬ 
dence that the mechanisms responsible for the X-ray/VHE emis¬ 
sion during non-flaring-activity states might not differ substan¬ 
tially from the ones responsible for the emission during flaring- 
activity states. In particular, the positive X-ray/VHE correla¬ 
tion observed during flaring activity in Mrk 421 and many other 
blazars has been interpreted by many authors as evidence for lep- 
tonic scenarios, and hence against the hadronic scenarios where 
the X-ray emission and the y-ray emission are produced by dif¬ 
ferent particle populations and processes. However, with a fine 
tuning of the parameters, hadronic models are also able to ex¬ 
plain single flaring events with an X-ray/VHE correlation with 
time lag zero (Mastichiadis et al.|2013|l. 


The observations presented here confirm that the relation be¬ 
tween X-ray and VHE bands also exists when the source is not 
flaring. That is, such a relation is persistent over long timescales 
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of at least several months, and does not occur only on single 
flaring events, and that is much more difficult to explain with 
hadronic scenarios. Therefore, these observations provide strong 
evidence supporting leptonic scenarios as responsible for the 
dominant X-ray/VHE emission from Mrk 421. 

Another result that is worth discussing is the overall anti¬ 
correlation between the X-ray and optical bands. This anti¬ 
correlation spreads over a large range of time lags and hence 
it is fundamentally different from the one obtained for the X- 
ray and VHE bands. As we show in Section]^ the origin of this 
overall anti-correlation is the long-term trends of the optical/UV 
and X-ray activity: while the former increases over time dur¬ 
ing the observing campaign, the latter decreases. The temporal 
evolution of the optical and X-ray/VHE bands is complex, and a 
dedicated correlation analysis over many years will be necessary 
in order to properly characterize it. 

Eor the 2009 multi-wavelength campaign, we do observe an 
anti-correlation in the (long-term) temporal evolution between 
optical/UV and X-rays, and hence it is worth discussing possible 
theoretical scenarios that might lead to this situation. The first 
scenario is that the optical/UV and the X-ray/VHE emissions 
are dominated by the emission from two distinct and uncon¬ 
nected regions with different temporal evolutions of their respec¬ 
tive particle populations. In such case, the optical/UV vs. X-ray 
anti-correlation observed in the 2009 multi-instrument campaign 
would have occurred by chance, and hence we would also expect 
to see multi-month time intervals with a positive correlation, or 
no correlation. A second scenario would be a two-component 
(high- and low-energy) particle population, in which the low- 
and high-energy particles have a different but related long-term 
temporal evolution. A change in the magnetic field intensity 
while keeping the acceleration timescale constant could lead 
to the observed optical/UV-X-ray (long-term) anti-correlation. 
An increase in the magnetization would produce a higher syn¬ 
chrotron emission with a decrease in the energy of the electrons 
(due to a stronger cooling), which effectively would lead to a 
higher optical emission with a lower X-ray emission. On the 
other hand, a lower magnetization would lead to a decrease in 
the total emitted synchrotron flux, but a higher maximum elec¬ 
tron energy, which effectively could lower the optical flux and 
increase the X-ray flux. In practice, such a scenario would be 
somewhat similar to the blazar sequence (Ghisellini et al.|1998 


with the difference that in the latter scenario the different coo - 
ings would relate to different sources instead of different states 
of the same source. A third scenario could be a global long¬ 
term change in the efficiency of the acceleration mechanism 
that produces the electron energy distribution. Such a change in 
the global efficiency could shift the entire synchrotron bump to 
higher/lower energies. Eor instance, if the acceleration mecha¬ 
nism becomes more efficient to get electrons up to the highest 
energies at the expense of keeping a lower number of low-energy 
electrons (i.e., the index of the electron population gets harder), 
the emission at the rising segment of the synchrotron bump (op¬ 
tical) would decrease, while that on the decreasing segment of 
the synchrotron bump (X-rays) would increase. 

In this study we did not see any correlation between the ra¬ 
dio fluxes and those at higher frequencies. However, since the 
measured radio emission is expected to have large contributions 
from regions farther away in the jet, it is not surprising to see 
a non-correlation between radio and optical. X-ray and/or y-ray 
energies on timescales of days to weeks. We note, however, that 
such correlation might be apparent during large flares when the 
radio emission might be strongly dominated by the same region 
responsible for the overall broadband emission. 


7. Conclusions 


We studied the broadband evolution of the SED of Mrk 421 
through a 4.5 month long multi-instrument observing campaign 
in 2009, when the source was in its typical (non-flaring) activ¬ 
ity state, with a VHE flux of about half that of the Crab Nebula. 
Even though the source did not show flaring activity, we could 
measure significant variability in the energy bands where the 
emitted power is largest: optical. X-ray, y rays and VHE. The 
highest variability occurred in the X-ray band, which, within 
the standard one-zone SSC scenario, indicates that the high- 
energy electrons are more variable than the low-energy elec¬ 
trons. We also observed a harder-when-brighter behavior in the 
X-ray spectra, and found a positive correlation between the X- 
ray and VHE bands. In the literature one can find many works re¬ 
porting a positive correlation between the X-ray and VHE fluxes 
(e.g. Fossati et al. 2008 and references therein) and spectral 
shape changes with the X-ray flux (e.g. |Tramacere et al. [2009) 1, 
but only when Mrk421 was showing VHE flaring activity (i.e., 
VHE flux above the flux of the Crab Nebula). This is the first 
time that such characteristics are reported for non-flaring activ¬ 
ity and suggests that the processess occurring during the flaring 
activity also occur when the source is in a non-flaring (low) state. 
In particular, this is a strong argument in favor of leptonic scenar¬ 
ios dominating the broadband emission of Mrk 421 during non¬ 
flaring activity, since such a temporally extended X-ray/VHE 
correlation cannot be explained within the standard hadronic 
scenarios. Moreover, a negative correlation in the (long-term) 
temporal evolution of the optical/UV and X-ray bands was also 
observed. Such a trend could be produced in a region with a par¬ 
ticle population where the low- and high-energy particles evolve 
differently but in a related way, which could be produced by 
a change in the magnetization of the region while keeping the 
acceleration timescales constant, or by a global change in the 
efficiency of the mechanism accelerating the electrons. In any 
case, even though statistically significant for the 2009 multi¬ 
instrument campaign, the current dataset does not allow us to 
exclude that this optical/X-ray anti-correlation was observed by 
chance, and hence that the optical and the X-ray bands are pro¬ 
duced by distinct and unrelated particle populations that evolve 
separately. Further multi-instrument observations extending over 
many years will help to address this question. 
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Fig.A.l. Measured fractional variability Fyar (black) compared 
with a delete-tf jackknife estimate (red) with d = W, for all 
light curves with a minimum of 15 data points. 
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Appendix A: Reliability of Fyar 

Almost all light curves of the campaign are unevenly sampled. 
The sampling is different for each instrument depending on ob¬ 
servation schedule, weather and technical issues. There are of¬ 
ten gaps of different lengths in the light curves and each light 
curve has a different number of data points ranging from a few 
up to a few hundred. In addition, some light curves are binned 
into bins of several days because of the limited sensitivity of the 
corresponding instruments. Therefore we have to assess whether 
there is an error introduced to by the uneven sampling and 
the binning and how the Fy^ values can be compared. In addition 
we need to know the minimum number of flux measurements per 
light curve which are needed to obtain a reliable, unbiased Fvar 
value. 

To address these questions, we hrst made a delete-^f jack¬ 
knife analysis, i.e., from each light curve containing N data 
points, we randomly removed d data points (I < d < N - 2) 
and calculated Tvar for this reduced sample ofN-d data points. 
Applying a delete-^f jackknife analysis on a time series is for¬ 
mally not correct, as the data are correlated in time (blazar light 
curves usually show a red-noise behavior). However, our pur¬ 
pose is to create datasets that have statistical properties identi¬ 
cal to the real data to demonstrate the impact of gaps and un¬ 
even sampling on Fy^^. Figure [AT| shows Fyar for the jackknife 
datasets with d - '^fN in comparison with Fyar for the original 
light curves. The Fyar values do not change significantly. The er¬ 
ror bars are larger because of the reduced number of flux values 
in the jackknife samples. The result does not depend on the par¬ 
ticular choice of d, as long as there are sufficient datapoints in the 
jackknife-samples remaining. Some of the light curves are (al¬ 
most) regularly sampled, namely Termi-LAT, RXTE/PCA and 
RXTE/ASM. These are good examples that demonstrate that ir¬ 
regular sampling does not introduce a bias to the Fyar measure¬ 
ment. 

We also varied d between 1 and N-2. Eigure [A^ shows Fy^^ 
vs. {N - d)IN for selected light curves. As long as A -is larger 
than ~ 5, the measured Avar is approximately constant with vary¬ 
ing d and in agreement with Fy^y of the original light curve. 
Likewise, the error bars do not change significantly. Strong de¬ 
viations from Avar of the original light curve occur only, if at 
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Fig.A.2. Fractional variability F^ar (red) as a function of {N - d)IN for the jackknife-t/ samples of selected representative light 
curves. The measured F^ar of the original light curve and its error are shown as red horizontal lines. The F^,ar of the jackknife- 
samples is constant and agrees with the original F^ar within its errors for all but the largest d (i.e., the smallest jackknife-samples). 



Fig.A.3. Measured fractional variability F^,ar (black) compared 
with a moving-block jackknife estimate (red), using a blocksize 
of m - ifN, for all light curves with a minimum of 15 data 
points. 


all, when N - d < 10. No significant deviations are observed 
at N - d > 20 for any of the light curves. Thus we conclude 
that our Fyar measurement is robust for all light curves with 20 
or more flux data points. In our multi-wavelength sample most 
light curves in the optical, UV, X-rays, HE y rays and VHE have 
more than 20 flux data points. In the radio and near-infrared, 
most light curves do not. 

Fvar is reliable for all but the smallest samples. If after re¬ 
moval of d data points the remaining sample is smaller than 
about 10, then Fvai might be over- or underestimated. 

We also made a moving-block jackknife test, i.e., we re¬ 
moved blocks of m consecutive measurements. This test is still 
formally correct when the data are slightly correlated in time, but 
the drawback is that the number of jackknife samples is much 
smaller than in case of the delete-c/ jackknife test. Likewise, it 
only shows the i nfluen ce of gaps, not the influence of a random 
sampling. Eigure [a!^ shows the test for m - Va. As in the case 
of the delete-t/ jackknife test, the Fyar values do not change sig¬ 
nificantly and the uncertainties are larger. 

These conclusions, however, are only valid because the 
dataset does not have any strong flare and hence it is unaffected 
by removing points or blocks randomly. Occasional strong flares 
therefore should be removed before doing an Fyar analysis of 
light curves that are otherwise in a typical or low state. 
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